Phase Separation in Electronic Models for Manganites 
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The Kondo lattice Hamiltonian with ferromagnetic Hund's coupling as a model for manganites is 
investigated. The classical limit for the spin of the (localized) t2 g electrons is analyzed on lattices 
of dimension 1,2,3 and oo using several numerical methods. The phase diagram at low temperature 
is presented. A regime is identified where phase separation occurs between hole undoped antiferro- 
magnetic and hole-rich ferromagnetic regions. Experimental consequences of this novel regime are 
discussed. Regions of incommensurate spin correlations have also been found. Estimations of the 
critical temperature in 3D are compatible with experiments. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 
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The phenomenon of colossal magnetoresistEtn.ee in 
metallic oxides such as Lai_ x Sr x Mn03 has recently at- 
tracted considerable attention fm due to its potential 
technological applications. A variety of experiments have 
revealed that oxide manganites have a rich phase di- 
agram Q with regions of antiferromagnetic (AF) and 
ferromagnetic (FM) order, as well as charge ordering, 
and a peculiar insulating state above the FM critical 
temperature, Tf M . Recently, layered manganite com- 
pounds La1.2Sr1.8Mn.2O7 have also been synthesized ||] 
with properties similar to those of their 3D counterparts. 

The appearance of ferromagnetism at low tempera- 
tures can be explained using the Double Exchange (DE) 
mechanism HJq]. However, the DE model is incomplete 
to describe the entire phase diagram observed experimen- 
tally. For instance, the coupling of the electrons with the 
lattice may be crucial to account for the insulating prop- 
erties above Tf M ||. The presence of a Berry phase 
in the large Hund coupling limit also challenges predic- 
tions from the DE model M. In this paper we remark 
that another phenomena occurring in manganites which 
is not included in the DE description, namely the charge 
ordering effect, may be contained in a more fundamen- 
tal Kondo model where the t^g (localized) electrons are 
ferromagnetically (Hund) coupled with the e g (mobile) 
electrons. More precisely, here we report the presence of 
phase separation (PS) between hole undoped antiferro- 
magnetic and hole-rich ferromagnetic regions in the low 
temperature phase diagram of the FM Kondo model. 
Upon the inclusion of long-range Coulombic repulsion, 
charge ordering in the form of non-trivial extended struc- 
tures may be stabilized. A similar phenomenon has been 
extensively discussed before in the context of the high-Tc 
superconductors, although without the presence of ferro- 
magnetic domains [^,^|. In the cuprates neutron scat- 
tering experiments have shown indications of stripe for- 
mation |l(| , as predicted by theoretical calculations after 
the inclusion of Coulomb interactions |ll| . The analysis 
reported in this paper suggests that phenomena as rich 
as observed in the cuprates may exist in the mangan- 



ites, and hopefully the present results will induce further 
theoretical and experimental work in this context. 
The FM Kondo Hamiltonian is defined as 



H 



-t 



E 



(ct 

\ 1(7 



h.c. 



ia/3 



0~a/3Ci/3 ■ Si, (1) 



where c\„ are destruction operators for one species of e g - 
fermions at site i with spin tr, and Si is the total spin 
of the ti g electrons, assumed localized. The first term 
is the e g electron transfer between nearest-neighbor Mn- 
ions, Jh > is the Hund coupling, the number of sites is 
L, and the rest of the notation is standard. The density is 
adjusted using a chemical potential /1. In this paper the 
spin Si will be considered classical (with |Si| = 1) rather 
than quantum mechanical, unless otherwise stated. Phe- 
nomenologically Jh 3> t, but here Ju/t was considered 
an arbitrary parameter, i.e. both large and small values 
for Jn/t were studied. Although models beyond Eq.(l) 
may be needed to fully understand the manganites, it is 
important to study the properties of simple Hamiltonians 
to clarify if part of the experimental rich phase diagram 
can be accounted for using purely electronic models. 

To study Eq.(l) in the t^g spin classical limit a Monte 
Carlo technique was used: first, the trace over the 
fcrmionic degrees of freedom in the partition function was 
carried out exactly diagonalizing the 2L x 2L Hamiltonian 
of electrons in the background of the spins {Si}, using li- 
brary subroutines. Here the fermionic trace is a positive 
function of the classical spins and, thus, the resulting in- 
tegration over the two angles per site parametrizing the 
Si variables can be performed with a standard Monte 
Carlo algorithm without "sign problems" . In addition, 
part of the calculations were also performed with the 
Dynamical Mean-Field approximation (D=oo) |l2| , the 
Density-Matrix Renormalization Group (DMRG), and 
the Lanczos method. Special care must be taken with 
the boundary conditions (BC). Closed shell BC or open 
BC are needed to stabilize a ferromagnetic spin arrange- 
ment. If other BC are used the spin correlations at short 
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distances are still strongly FM (if working at couplings 
where ferromagnetism is stable), but not at large dis- 
tances where they become negative. This well-known ef- 
fect was observed before in this context [|l3],[l4| and it 
does not present a problem in the analysis shown below. 



15.0 



10.0 



5.0 



0.0 



15.0 



10.0 



5.0 



0.0 









D = 1 




(a) 






i 1 FM 








PS \. 















D = 2 




(b) 




FM 


ps\ 









0.02 




0.00 



1.0 0.8 0.6 0.4 0.2 0.0 1.00 0.95 0.90 0.85 
<n> <n> 

FIG. 1. Phase diagram of the FM Kondo model reported 
here. FM, IC, PM and PS denote regimes with strong FM 
correlations, incommensurate correlations, paramagnetic cor- 
relations, and with phase separation between undoped AF 
and hole-rich FM regions, respectively, (a) was obtained with 
MC simulations at T — t/75 using chains with L=20, 30 and 
40 sites. The full circles correspond to classical spins. The 
stars and dashed lines are DMRG results obtained with t2 g 
spins S=3/2 on chains with up to L=16 sites keeping 48 states, 
and using (2/3) Jh as coupling in the Hund term; (b) are MC 
results for classical spins at T = t/50 using 4 x 4, 6 x 6 and 
8x8 clusters. In the shaded region a crossover from PS to IC 
was observed but the actual boundary position is difficult to 
find accurately; (c) corresponds to results in the D=oo limit 
at a fixed coupling Jh/W — 4.0 and as a function of T and 
(n) (W is defined in the text). The "PS (AF-PM)" region 
denotes PS between undoped AF and hole-rich PM regions. 



The results reported here can be summarized in the 
phase diagrams obtained for several dimensions D shown 
in Fig.l. In both ID and 2D and at low temperatures 
clear indications of (i) strong ferromagnetism, (ii) incom- 
mensurate (IC) correlations, and (iii) a novel regime of 
phase separation were identified. For D— 1 and 2, finite 
size effects were found to be small for the lattice sizes used 
in this study, although the PS-IC boundary in 2D was 



difficult to identify accurately and, thus, only a crossover 
is indicated. Results arc also available in small 3D clus- 
ters and qualitatively they agree with those in Figs.la-b. 
In ID we also obtained results with quantum t^g spins 
S=3/2. In this case the PS regime was studied calculat- 
ing the compressibility with the DMRG technique, and 
the FM regime monitoring the ground state spin quan- 
tum number with the Lanczos method. The agreement 
with results in the classical limit is good, specially re- 
garding the ferromagnetic line. PS for S=3/2 is between 
hole undoped AF and hole-rich non-fully-saturated FM 
regions. Similar results were obtained in the case S=l/2. 
In Fig.lc results in D=co are shown at Jh/W — 4, where 
W is the half-width of the semicircular density of states 
D{e) = (2/7iTT)^l - (e/W) 2 for the e g electrons. In 
excellent agreement with the predictions for D=l and 
2, at low temperatures phase separation between un- 
doped AF and hole-rich FM regions was observed. For 
(ri) < 0.88 and T/W -C 1 the ground state becomes 
ferromagnetic. The quantitative similarities between the 
results obtained for D=l,2 and oo led us to believe that 
the conclusions of this paper are also valid for 3D man- 
ganite systems. Below the details justifying the proposed 
phase diagrams shown in Fig.l are provided. 
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FIG. 2. Spin-spin correlations of the classical spins at zero 
momentum S(q = 0) vs temperature T (in units of i). Results 
for several lattice sizes are shown on (a) chains and (b) 2D 
clusters. The technique used is the Monte Carlo algorithm 
described in the text. The density and coupling are shown. 
In (a) closed shell BC are used i.e. periodic BC for L = 10 
and 30 and antiperiodic BC for L = 20 and 40. In (b) open 
BC are used. 

The boundaries of the FM region of the phase di- 
agram were found evaluating the spin-spin correla- 
tion between the classical spins defined as 5(q) = 
( 1 /i)Ej, m e l(j " m) q (S j • S m >. Fig.2 shows 5(q) at zero 
momentum vs. T/t for typical examples in ID and 2D. 
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The rapid increase of the spin correlations as T is reduced 
and as the lattice size grows clearly points towards the 
existence of ferromagnetism in the system . Repeat- 
ing this procedure for a variety of couplings and densi- 
ties, the robust region of FM shown in Fig.l was deter- 
mined, in qualitative agreement with previous studies at 
D = oo |Q and in the large Jh /t limit on ID chains [ p^[ . 

In the small Jjj/t region incommensurate correlations 
were observed. IC effects were predicted using a Hartree- 
Fock approximation (ll| as an alternative to canted 
FM [§. The IC region was observed in our studies mon- 
itoring 5(q). Both in ID and 2D there is one dominant 
peak which moves away from the AF location at (n) = 1 
towards zero momentum as (n) decreases. In the 2D 
clusters the peak moves from (tt, tt) towards (7r, 0) and 
(0, 7r), rather than along the main diagonal. Our results 
confirm a tendency towards the development of IC cor- 
relations in the FM Kondo model However, our 
computational study predicts IC correlations only in the 
small and intermediate Jn/t regime i.e. in a region of 
parameter space not realized in the manganites. Details 
will be provided elsewhere. 
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FIG. 3. Electronic density (n) vs chemical potential \x in 
(a) D=l, (b) D=2, and (c) D = oo clusters. The temperatures 
are indicated. The coupling is Jn/t = 8.0 in (a) and (b) and 
Jh /W = 4.0 in (c). The meaning of W is explained in the 
text. PBC were used both in D=l and 2. The discontinuities 
shown in the figures are indicative of PS. In (a) the inset 
contains the spin-spin correlation in real space at densities 
1.0 and 0.72 showing that indeed PS occurs between AF and 
FM regions. 



The main result of the paper is contained in Fig. 3 
where the computational evidence for the existence of 
phase separation in dimensions 1, 2, and oo is given [fill. 
The presence of a discontinuity in (n) vs u, shows that 
some electronic densities can not be stabilized. If the 
system is nominally prepared with such density it will 
spontaneously separate into two regions with the densi- 
ties corresponding to the extremes of the discontinuities 
of Fig. 3 By analyzing these extremes the proper- 

ties of the two domains can be studied. At least in the 
classical limit for the ti g spins, one region is undoped 
((n) = 1) with strong AF correlations, while the other 
contains all the holes and the spin-spin correlations be- 
tween the classical spins indicate the presence of strong 
FM correlations [^] (for a justification of this statement 
see the inset of Fig. 3a. The results are similar in D=2 
and infinite). This is natural since holes optimize their 
kinetic energy in a fully aligned background. On the 
other hand, at (n) — 1 the DE mechanism is not oper- 
ative: if the electrons fully align their spins they simply 
cannot move in the conduction band due to the Pauli 
principle. Then, energetically it is better to form an an- 
tiferromagnetic pattern (similar intuitive and numerical 
conclusions were reached for the case of S=3/2 and 1/2 
localized spins). As Jh grows, the jump in Fig. 3 is re- 
duced and it tends to disappear in the Jh = oo limit 
& 

Experimentally, PS may be detectable using neutron 
diffraction techniques if the two coexisting phases have 
different lattice parameters as in La2Cu04+5, a Cu-oxide 
with hole-rich and hole-poor regions [Q. NMR and 
NQR spectra, as well as magnetic susceptibility mea- 
surements, can also be used to detect PS p3|^4| . Note 
also that in the PS regime 5(q) presents a two peak 
structure, one located at the AF position and the other 
at zero momentum. Since this also occurs in a canted 
ferromagnetic state care must be taken in the analy- 
sis of the experimental data. In particular recent ex- 
perimental results by Kawano et al. |25f| are in qualita- 
tive agreement with Fig.lc since these authors observed 
a reentrant structural phase transition accompanied by 
"canted ferromagnetism" below T C FM , at 0.10 < x < 0.17 
in Lai_ x Sr x Mn03. Also the polaron-ordered phase re- 
ported by Yamada et al. p6| ] can be reanalyzed in terms 
of the present results since it is known that the AF 
phase in 3D manganites is orthorhombic while the FM is 
pseudo-cubic. The formation of a lattice superstructure 
may stabilize the magnetic tendency to phase separate 
and minimize lattice distortions. 

Another alternative is that Coulombic forces and low 
Sr mobility prevent the macroscopic accumulation of 
charge intrinsic of a PS regime. Thus, extended hole- 
rich domains in the form of stripes (as in cuprates) [jll| , 
or some other arrangement, could be formed. Although 
these details certainly deserve further work, the results 
in this paper are enough to show that tendencies similar 
to those found in models for the 2D cuprates, notably 
the t — J model, may also be operative in manganites. 
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The main difference is that instead of separation between 
hole-poor AF and hole-rich superconducting regions, as 
in La2Cu04+5, here the observed PS is between undoped 
AF and hole-rich ferromagnetic regimes. Tendencies to 
pairing in Cu-oxides are replaced by tendencies to form 
ferromagnetic ground states in Mn-oxides (ll|] . 

Although the phase diagrams of Fig.l have PS close 
to half-filling, actually this phenomenon should also oc- 
cur at high hole doping (n) ~ if an extra direct AF 
exchange interaction between the localized ti g spins is 
included. This coupling may be originated in a small 
hopping amplitude for the ti g electrons. At (n) = 
model Eq.(l) supplemented by a Heisenberg coupling 
J'/t among the localized spins will certainly produce an 
AF phase, as in experiments, which upon electron doping 
will induce a competition between AF (with no e g elec- 
trons) and FM electron-rich regions, similarly as in the 
half-filled limit but replacing holes by electrons. Previ- 
ous studies in ID support these claims Jll|. Thus, PS or 
charge ordering could exist in manganites both at large 
and small fermionic densities. 
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FIG. 4. (a) Real space spin-spin correlations among the 
classical spins at distance \i — j\ obtained with the MC tech- 
nique at (n) =0.5, Jh = oo, on a 6 3 cluster, and parametric 
with temperature T/t. The inset shows S(q) at zero momen- 
tum and several densities vs T/t; (b) Bounds on the critical 
temperature as estimated from data as shown in (a) as a func- 
tion of (n). 

For completeness, upper bounds on the critical temper- 



ature T C FM in 3D systems are also provided. Using Monte 
Carlo simulations in principle it is possible to calculate 
Tf M accurately. However, the characteristics of the al- 
gorithm used here prevented us from studying clusters 
larger than 6 3 even at Jh = oo. In spite of this limita- 
tion, monitoring the spin-spin correlations in real space 
allows us to judge at what temperature T* the correla- 
tion length reaches the boundary of the 6 3 cluster. Since 
the bulk T FM is smaller than T*, this allows us to es- 
tablish upper bounds for the critical temperature based 
on finite size cluster calculations. Fig. 4a shows the spin- 
spin correlations at several temperatures in the Jh = oo 
limit and at (n) = 0.5. When T ~ O.lt robust correla- 
tions reach the boundary, while for T > 0.12i the corre- 
lation is short-ranged. Thus, at this density we estimate 
that Tf M < 0.12i. Results for several densities in this 
limit are shown in Fig. 4b. For realistic densities, such 
as (n) ~ 0.7, our results are about a factor 1.7 smaller 
than predicted by high temperature expansions (this 
small discrepancy may be caused by the use of S=l/2 
localized spins in Ref. |p7| ). Nevertheless the order of 
magnitude of both calculations is similar. Monitoring 
the rapid growth of the zero momentum spin correla- 
tions (inset of Fig. 4a) provides similar bounds. T C FA/ is 
the highest at (n) — 0.5, if Jh — oo. Since results for 
the e g electrons bandwidth range from BW ~ 1 eV p3] 
to BW ~ 4eV producing a hopping t = BIT/12 be- 
tween 0.08 and 0.33 eV, then our estimation for the crit- 
ical temperature ranges roughly between Tf M ~ 100 K 
and 400 K. This is within the range observed exper- 
imentally, in agreement with other results |l^,^], and 
in disagreement with previous estimations using classical 
spins that predicted a much higher critical temperature 
Tf M - 2000X @. Then, purely electronic models can 
account for Tj without resorting to electron-phonon 
couplings. 

Summarizing, here the main conclusion of the paper 
is that the phenomenon of phase separation occurs in 
realistic models for manganites. Then, experimentalists 
shoud consider this potential regime in the interpreta- 
tion of their results i.e. tendencies to phase separate 
could be as strong as those towards ferromagnetism in 
the real materials. Hopefully the present paper will trig- 
ger discussion on how phase separation will reveal itself 
in manganites from the experimental point of view. It 
is anticipated that the inclusion of long-range Coulombic 
interactions may stabilize complex charge ordered struc- 
tures, that could be observed experimentally using a vari- 
ety of techniques. Neutron scattering experiments should 
detect these charge fluctuations above the critical tem- 
perature, similarly as it occurs in the cuprates. From a 
more basic perspective, our results suggest that the ex- 
istence of phase separation in the context of electronic 
models for transition metal oxides is a more general phe- 
nomenon than previously anticipated. 

S. Y. is supported by the Japanese Society for the Pro- 
motion of Science. J. H. is supported by the Florida State 
grant E&G 502401002. He thanks E. Miranda for useful 



4 



discussions. A. L. M. acknowledges the financial support 
of the Conselho Nacional de Desenvolvimento Cientffico 
e Tecnologico (CNPq-Brazil) A. M. and E. D. are sup- 
ported by the NSF grant DMR-9520776. 



[1] S. Jin, et al., Science 264, 413 (1994). 

[2] P. E. Schiffer, A. P. Ramirez, W. Bao, and S-W. Cheong, 

Phys. Rev. Lett. 75, 3336 (1995); C. H. Chen and S-W. 

Cheong, Phys. Rev. Lett. 76, 4042 (1996). 
[3] Y. Moritomo, et al., Nature 380, 141 (1996). 
[4] C. Zener, Phys. Rev. 82, 403 (1951); P. W. Anderson 

and H. Hasegawa, Phys. Rev. 100, 675 (1955). 
[5] P. G. de Gennes, Phys. Rev. 118, 141 (1960). 
[6] A. J. Millis, P. B. Littlewood, and B. I. Shraiman, Phys. 

Rev. Lett. 74, 5144 (1995); H. Roder, J. Zang and A. R. 

Bishop, Phys. Rev. Lett. 76, 1356 (1996). 
[7] E. Muller-Hartmann and E. Dagotto, Phys. Rev. B 54, 

R6819 (1996). 

[8] V. J. Emery, S. A. Kivelson, and H. Q. Lin, Phys. Rev. 
Lett. 64, 475 (1990). Note that here PS between AF hole- 
poor and FM hole-rich regions was also mentioned. How- 
ever, in the t — J model this may occur at unrealistic 
values of the couplings, while in the FM Kondo model 
it is present in the realistic regime. We thank S. Kivel- 
son for this comment. See also V. J. Emery, and S. A. 
Kivelson, Physica C 209, 597 (1993). 
[9] E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

[10] J. M. Tranquada et al., Nature 375, 561 (1995). 

[11] U. Low et al, Phys. Rev. Lett. 72, 1918 (1994); S. Haas 
et al., Phys. Rev. B 51, 5989 (1995). 

[12] N. Furukawa, J. Phys. Soc. Jpn. 63, 3214 (1994); ibid, 
2754 (1995). 

[13] J. Riera, K. Hallberg and E. Dagotto, Phys. Rev. Lett. 
79, 713 (1997). 

[14] K. Kubo, J. Phys. Soc. Jpn. 51, 782 (1982); J. Zang, 
H. Roder, A. R. Bishop, and S. A. Trugman, J. Phys.: 
Condens. Matter 9, L157 (1997); T. A. Kaplan and S. D. 
Mahanti, J. Phys.: Condens. Matter 9, L291 (1997). 

[15] In both ID and 2D the Mermin- Wagner theorem forbids 
a nonzero T^ M in the model. However, the correlation 
lengths can be very large even at finite T/t. 

[16] J. Inoue and S. Maekawa, Phys.Rev.Lett.74, 3407 (1995). 

[17] Our results for the IC correlations in 2D are in qualitative 
agreement with M. Hamada and H. Shimahara, Phys. 
Rev. B 51, 3027 (1995). 

[18] Up to 10 6 sweeps were needed to obtain accurate results 
in Fig.3.a-b, maintaining the acceptance rate close to 0.5. 

[19] Actually, by fine tuning jj, to be very close to its "critical" 
value, a two peak structure was observed in (n) as in a 
first order phase transition. 

[20] It was observed that the FM correlations at the bound- 
ary with PS in Fig.l are weaker than at the center of 
the FM phase, probably caused by a small T — order 
parameter. 



[21] PS in manganites was also conjectured in Ref. jl3j based 

on ID studies of the H — J" version of Eq.(l). 
[22] P. G. Radaelli, et al., Phys. Rev. B 49, 6239 (1994) and 

references therein. 
[23] P. C. Hammel, et al., Phys. Rev. Lett. 71, 440 (1993). 
[24] It is interesting to observe that experimental studies 

of Sr2- y La y Mn04 have reported evidence for electronic 

phase separation and charge ordering at small y. See W. 

Bao, et al., preprint. 
[25] H. Kawano et al., Phys. Rev. B 53, R14709 (1996); Phys. 

Rev. B 53, 2202 (1996). 
[26] Y. Yamada et al., Phys. Rev. Lett. 77, 904 (1996). 
[27] H. Roder, R. R. P. Singh, and J. Zang, Phys. Rev. B 56, 

5084 (1997). 

[28] Y. Moritomo, A. Asamitsu and Y. Tokura, Phys. Rev. B 

51, 16491 (1995). 
[29] D. D. Sarma et al., Phys. Rev. B 53, 6873 (1996). 



5 



